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Abstract 

We compute the relation between the quark mass defined in the min- 
imal modified MS scheme and the mass defined in the "Regularization 
Invariant" scheme (RI), up to the NNLO order. Lhe RI scheme is conve- 
niently adopted in lattice QCD calculations of the quark mass, since the 
relevant renormalization constants in this scheme can be evaluated in a 
non-perturbative way. The NNLO contribution to the conversion factor 
between the quark mass in the two schemes is found to be large, typically 
of the same order of the NLO correction at a scale [i ~ 2 GeV. We also 
give the NNLO relation between the quark mass in the RI scheme and the 
renormalization group-invariant mass. 



1 Introduction 



The values of quark masses are of great importance in the phenomenology of 
the Standard Model and beyond. For instance, bottom and charm quark masses 
enter significantly in the theoretical expressions of inclusive decay rates of heavy 
mesons, while the strange quark mass plays a crucial role in the evaluation of the 
K — > 7T7T, AI = 1/2, decay amplitude and of the CP- violation parameter e'/e. 
In the Standard Model, quark and lepton masses are fundamental parameters. 
Quark masses, however, cannot be measured directly, since quarks do not appear 
as physical states and are confined into hadrons. Thus, for quarks the kinematical 
concept of on-shell mass is meaningless. The values of the quark masses then 
depend on precisely their definitions, which can can be given in terms of a short 
distance mass in some renormalization scheme. 

A popular definition of the quark mass is the MS mass. For each flavour of 
quarks, this quantity is a well defined, short distance running coupling, provided 
the relevant scale /i is chosen inside the perturbative region of QCD. An alterna- 
tive, but equally convenient definition of the quark mass is the so called invariant 
quark mass. This has the major advantage of being both a scale and a scheme 
independent quantity. At present, the connection between the invariant mass and 
the MS mass is known up to four-loop order in QCD [1J. 

For heavy quarks only, the pole mass has also been considered and its relation 
to the MS mass has been computed up to two-loop order [Q. Since, however, for 
confined quarks there should not be any pole in the full propagator, the definition 
of the pole mass makes sense only in perturbation theory. On the other hand, 
the pole mass is affected by a renormalon ambiguity |3], |3J], which prevent the 
possibility of precisely define its perturbative expansion. For this reason, and 
also because the pole mass can be only defined in the heavy quark short 
distance definition of the quark mass, as the MS or the invariant quark masses, 
should be always preferred. 

Despite the great theoretical effort which has been spent in the last few years 
to determine the values of quark masses, these quantities still remain among the 
least well known parameters in the Standard Model. This is expecially true in the 
light quark sector, where the relative uncertainties on the up, down and strange 
quark masses are estimated to be as large as 50% 0. 

Lattice QCD is in principle able to predict the mass of any quark by fixing, 
to its experimental value, the mass of a hadron containing a quark with the same 
flavour ||. The quark mass that is directly determined in lattice simulations is the 
(short distance) bare lattice quark mass m(a), where a is the lattice spacing, the 
inverse of the UV cut-off. The conversion factor relating m(a) to a continuum 
renormalized mass, m(fi), only depends on typical scales of the order of \i ~ 
cr 1 ~ 2 — 4 GeV. In this respect lattice QCD is unique, since QCD sum rules 
calculations of quark masses have to work at much smaller scales, where higher- 
order corrections and/or non-perturbative effects |7j may be rather large. 
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It has been observed that perturbative calculations, within the lattice regular- 
ization, suffer from large higher-order corrections ||. This is due to the presence 
of tadpole-like diagrams, which are absent in continuum perturbation theory. 
Such corrections might then introduce a large uncertainty in the calculation of 
the renormalization constant relating the bare lattice quark mass m(a) to the 
renormalized mass m(n) in a continuum scheme. This constant is only known at 
one-loop order || [CO], and its NLO renormalization group-improved version has 



been derived in ref. 11 



The uncertainties related to the perturbative calculation of the lattice renor- 
malization constants can be completely avoided by using the non-perturbative 



renormalization program proposed in ref. ||12|| . Within this approach, the renor- 
malization constants of lattice bare correlation functions, or couplings and masses, 
are computed in a non-perturbative way in the so called "Regularization Invari- 
ant" (RI) scheme. When required, the connection between the renormalized 
quantities in the RI scheme and those defined in any other scheme, like the MS 
one, can be then computed in continuum perturbation theory. 

Recently, an extensive lattice calculations of light quark masses, with non- 
perturbatively determined renormalization constants in the RI scheme, has been 
performed in ref. [13[].[] The relation between the quark mass in the RI and MS 
schemes has been computed at one loop in ref. ||12|| . This result, combined with 



the two-loop expansion of the QCD /^-function and the quark mass anomalous 
dimension, has been then converted into a NLO renormalization group-improved 
determination of the quark masses. 

At present , both the /3-function and the mass anomalous dimension in the MS 
scheme have been computed up to four loops, in refs. |T3| and [III [T3| respectively. 



Thus, the two-loop conversion factor R m , relating the quark masses in the RI 
and MS schemes, is the only missing ingredient necessary to convert the non- 
perturbative lattice determinations of the quark masses to their MS counterparts 
at the NNLO. Such a calculation is the main result of this paper. The factor 
R m is given in eq. flT%|) for N c = 3 and in the Landau gauge, which is the gauge 
usually adopted in Rl-scheme lattice calculations. We find that the size of the 
NNLO corrections is quite large, of the order of 7% at a typical scale /i ~ 2 GeV, 
and comparable with the NLO contribution. 

The paper is organized as follows. In sec. |2|, we discuss the renormaliza- 
tion conditions for the quark mass, by referring in particular to the MS and RI 
schemes. The perturbative expressions for the relevant renormalization constants 
and the ratio i? m are given in sec. |3|. In sec. |] we discuss the evolution of the 
running quark mass, we introduce the invariant quark mass and express the result 
of our calculation as a relation between the invariant mass and the mass in the 
RI scheme. 



1 An explorative study can be found in ref. fjlfl. 
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2 Quark mass renormalization schemes 



The quark mass renormalization constant, Z m (/i), relates the bare quark mass 
m to the renormalized mass m (//), defined at a scale /i in a given scheme: 



(1) 



This definition provides the relation between the values of quark masses in dif- 
ferent renormalization schemes. By considering for definitiveness the MS and RI 
schemes, one finds: 

JUrt = =& = fc (2) 



m RI (/i ) Z MS 



Even though (/i) and (/x) are separately divergent in the infinite cut-off 
limit, the ratio R m is finite, being equal to the ratio of renormalized quark masses 
in different schemes. 

The quark mass renormalization can be conveniently defined through a renor- 
malization condition for the inverse bare quark propagator in momentum space, 
Sq 1 (p). In perturbation theory, Sq 1 (p) has the form: 



tS 1 (p)=^ 1 (p)-m ^ 2 (p) 



(3) 



where Si ;2 are (bare) scalar quantities which depend also on the cut-off and 
the bare coupling constant, mass and gauge parameter, a , m and £o- Once 
expressed in terms of renormalized parameters, a s , m and £, the renormalization 
of the quark propagator still requires the renormalization of the quark field itself. 
Thus, the renormalized propagator can be written as: 



iS 1 {p) = Z q j)Si(p)-Z m mS 2 (p) 



(4) 



where £^2 are the bare functions £1,2 expressed in terms of renormalized param- 
eters. 

Different renormalization prescriptions on £ lj2 define the expressions of the 
quark field and mass renormalization constants, Z q and Z m , in the correspond- 
ing schemes. In perturbation theory, these prescriptions can be conveniently 
expressed by giving the (finite) coefficients entering in the perturbative expan- 
sions of the renormalized quantities. Thus, for instance, for £ 12 , up to O (ccf), a 
renormalization scheme can be defined by: 



lim — Z„Tr 

m ^o48 q 

lim — Z q Z m Tr 



-0 12 



S 2 (p) 



1 + 



a* 



l r=-\if 



r 1 + 



(4tt) (4tt)' 



Si + . 



(5) 



a. 



p =- 



» 2 = l + W) 



r 2 + 



at 



(4tt)< 



S 2 + ... 
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where traces are taken over the color and spin indices. 

The popular MS scheme |16] amounts to require that, in dimensional regu- 
larization, the coefficients of the perturbative expansion of the renormalization 
constants only contain singular terms in 1/e, where D = 4 — 2e is the space-time 
dimension. Here we consider its standard modification, the MS scheme ||17|| . We 
choose to renormalize the gauge coupling constant a s and the gauge parameter £ 
always in this scheme, independently on the renormalization condition adopted 
for the quark field and mass. A notable feature of the MS scheme is that the 
renormalization constants of gauge-invariant operators and parameters, like the 
quark mass itself, are gauge independent. Consequently, the renormalized quark 
mass, in the MS scheme, is gauge independent. 

One of the advantage of the RI scheme is that the corresponding renormal- 
ization prescriptions can be stated in a non-perturbative way. In the case of S lj2 , 
for instance, one simply requires the left-hand sides of eqs. (|5|) to be equal to 
unity. In perturbation theory, this implies the vanishing of all the coefficients 
and Sj. The renormalization conditions in the RI scheme clearly depend on the 
choice of the gauge parameter £*, and different choices define in fact different RI 
renormalization schemes. Throughout this paper we will present results for the 
RI scheme in a generic covariant gauge, even though a common choice in this 
context is the Landau gauge, corresponding to £* = 0. Indeed, this choice is 
particularly convenient to implement in non-perturbative calculations of lattice 
QCD. 

In numerical lattice simulations, the quark field renormalization condition 
of eq. (|D is not easy to implement, because of the presence of the continuum 
derivative. A common practice is to define a different quark field renormalization 
through the condition: 



lim — Z'Tr 

12 9 



El (p) p 2 = _ /J 2 = 1 (6) 



In the Landau gauge, the difference between the two definitions appears at order 
O(o,g), and it must be taken into account at the NNLO accuracy. For this reason, 
we have also computed the ratio of the two relevant renormalization constants 
and the result will be given in the next section. 

In the calculation of S 12 an important check is provided by the vector chiral 
Ward identity, which relates these functions to the amputated Green functions of 
the vector current and scalar density respectively between external quark states, 
IV (p) and Ts(p)- In particular, in the chiral limit: 



fr) 

dPn 



s 2 (p) = r 5 ( P )j m=o (7) 



m=0 



We have verified that the above equations are indeed satisfied at the level of 
bare Green functions. By requiring the Ward identities to be satisfied also by 
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renormalized quantities, one obtains the relations: 

Zy = 1 , Z m = Zg 1 (8) 

among the renormalization constants in a generic scheme. In the MS— NDR 
scheme, the vector and axial Ward identities are unaffected by the minimal sub- 
traction procedure. Thus, in this scheme, eqs. (|8|) are automatically satisfied. In 
the RI scheme, in order to satisfy the Ward identities, the following renormaliza- 
tion prescriptions on rv and Ts must be consistently applied: 



lim -^Z RI 

m^o 48 q 



Tr 



lim — Z„ RI (Z RI V ^rT 

ti->o 12 q V 6 > 



7jJV„(p) 



(9) 



The above conditions can be also applied at finite values of the quark mass, 
provided a sufficiently large value of the renormalization scale is chosen, /i 2 ^> m 2 . 



3 Results for the renormalization constants 

In order to present the results of our calculation, we expand a generic renormal- 
ization constant Z as a series in the strong coupling constant: 

Z = 1 + ^z« + ^U< 2 ) + . . . (10) 

47T (47r) 2 v ' 

Subsequently, each coefficient is expanded in inverse powers of e: 

IV 

3=0 



By definition, all the coefficients Zq vanish in the MS scheme. In addition, by 
requiring that the ratio of renormalization constants in different schemes, namely 

Z RI /Z MS , is finite, one finds a set of useful identities between the zj s in the two 
schemes. These observations, combined together, imply the following relations: 

msV 1} _ f?m\u 



z M S y 2) = ( z RI)( 2 ) (12) 

z Msy 2) = (zM)f -(zM)f (zW) o w 
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s 3 s 4 

Figure 1: One- and two-loop diagrams for the quark propagator. 



Thus, the renormalization constant in the MS scheme can be completely derived 
from its counterpart in the RI scheme (but not vice- versa). One also finds that 

the perturbative expansion of the ratio Z^/Z^^ is simply given, up to two- loop 
order, by the finite coefficient of perturbative expansion of Z : 

Rm ~z^- 1+ W)y m )o + j^f{ Zm )o + - (13) 

m 

We have calculated the quark field and mass renormalization constants, Z q 
and Z m , in the MS and RI schemes. Adopting naive dimensional regularization, 
in a generic covariant gauge, we have computed the Feynman diagrams shown in 
fig. [L|. Since the corresponding amplitudes can be expanded up to linear terms 
in the bare quark mass, this only requires the evaluation of one- and two-loop p- 
integrals, which we have computed by using the method of integration by parts of 
ref. |y|. The Ward identities of eqs. (0) have been used to check the correctness 



of our calculation. For the quark field renormalization constant, in the RI scheme, 
we find: 



Z 



•RI\ (1) = (N* - 1) 

? J. 1 AT [ ^ ' 



(1) 


(iV 2 - 1) 


1 


2N C 


(1) 


(iV 2 - 1) 





AN C 


(2) 


(iV 2 - 1) 


2 


8iV c 2 


(2) 


(iV 2 - 1) 


1 


16iV c 2 



z « Jo " 4iV c 

zf 1 )^ = ^|^e(-e + 3iV 2 + 2eiV 2 ) (14) 

= AT^ (~ 3 - 2 ^ - 22N c - 8 ^ N c + + 4 ^j) 
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Z?)™ = (l-6e 2 -H5iV 6 2 - 76eiV c 2 -6e 2 iV e 2 

+48 £ AT 2 C 3 + 48 N?( 3 + 20 N c n f 

where ( is the Riemann zeta function and £3 = 1.20206. . .. The quark mass 
renormalization constant, in the same scheme, is given by: 

Z^ 1 ) = ^71^ (-3) 
Z% 1 ) ^TT^-i-S-H) 



(1) 


(iV 2 - 1) 


1 


2N C 


(1) 


(N? - 1) 





AN C 


(2) 


(N? - 1) 


2 


8 AT? 


(2) 


(iV 2 - 1) 


1 


48iV c 2 


(2) 


(iV 2 - 1) 





96iV c 2 



(^)f = ^^(-9 + 31^-4JV e n/) (15) 
Z ™ l T = (~ 135 - 54 e - 59 N 2 C + 54 g iV 2 + 20 N c n f ) 



Z^ 1 ) = [ n l Ar9 ' (-75 - 144 £ - 36 £ 2 - 2645 N 2 C - 108 £ iV 2 

-18 £ 2 iV 2 + 288 Cs + 576 N 2 C ( 3 + 356 N c n f 

The corresponding values of renormalization constants in the MS scheme can be 

obtained from eqs. (|T4D and (|TJ) by using eq. flT2"|). The result for is in 

agreement with the original calculation of ref. [19]. Notice that, as expected, 

Z^j~ turns out to be gauge- independent. 

We have also calculated the ratio A q between the quark field renormalization 
constants and Z' q , the latter being defined in eq. (^). The result, expanded 
in series of the strong coupling constant, is: 

A = 5^ = 1 + -^A (1) + a * A (2) + (16) 
Z> + (47r) A? + (4vr) 2 q + '" 1 ' 

where we find: 

A« = (17) 

Af = ^°~a (3 " + 22iV 2 + 14£iV 2 + 4£ 2 iV 2 - 4iV> , 

Notice that, in the Landau gauge (£ = 0), at a typical scale /i ~ 2 GeV, A q 
represents a tiny correction, smaller than 1%. 

Finally, we discuss the results for the ratio R m , which provides the relation 
between the quark mass in the MS and RI schemes at a fixed scale [i. This ratio 

is obtained by substituting the values of {z~~} and {z^~) lYl ec L- (HID- For 



convenience, we present here its numerical expression as obtained for N c = 3 in 
the Landau gauge: 

m 3 (4tt) V 9 3 s3 9 V (4tt) 2 

Eq. ( Jl8| ) is the main result of this paper. Numerically, one finds that the size of 
the NNLO contribution to R^ N , at a scale \i =2 GeV (rif = 4), is about 7%. 
This represents a significant correction, comparable to the NLO contribution. 



Quark mass evolution and the invariant quark 
mass 



In this section we discuss the NNLO relations between quark masses in different 
renormalization schemes and at different renormalization scales. This discussion 
also provides us with the opportunity to introduce the so called renormalization 
group-invariant quark mass, m, and its perturbative relation with the quark mass 
in the RI scheme. 

In a generic renormalization scheme, the renormalized quark mass obeys the 
following renormalization group equation: 



d 



lr, 



dfi 2 + 2 



P(a s ) 



_d_ 

da. 



+ 



f \ C 9 , 7m 



m(/i)=0 (19) 



Both the /3-function and the mass anomalous dimension in the MS scheme have 
been computed up to four loops in |14j] and Jl|, |15j respectively. In the NNLO ap- 
proximation we are considering here, we only need the corresponding expansions 
up to three loops. The QCD /3-function is given by: 



47T 



_d_ 

dfi 2 



Pi 
MS 



ll Ar 2 

— N c n f 

3 3 1 



(20) 



(iV c 2 -l), 



2857 



-N* + 



(N> - if 



54 

1415 Ar2 11 (N 2 C 

N 2 n f + — ' 

54 c } 18 



205 
136" 



N 



Urn 



79 



-ni H N c ni 

f 54 / 



The mass anomalous dimension in the MS scheme is: 



7K 



2Z- l a 2 



d 
dfi 2 



h V4vr 



8 



' m 



N c 



(i) 



iV„ 2 



1 



7 



T (2) 



N?-l 



3 203 Ar2 5 Ar 
— + iV 2 - -N c n f 

4 12 c 3 ' 



+71/ 



129 129 Ar2 11413 Ar4 

+ n: 



108 



23 



-N r . 



1177 



2"" 54 " c 



iV 3 - 12iV c C 3 - 12iV c 3 C; 



The evolution of the renormalized gauge parameters is given by: 



4 



(0) 



N r . /13 



(i) /ttsV 



i=l 



47T/ 



(21) 



(22) 



For completeness, we also present the coefficients of the perturbative expan- 
sion of the quark field anomalous dimension, up to the NLO. Using the results of 
eq. ([14]), we derive, in the MS scheme and in a generic covariant gauge: 



3 + 22N'Z + 8£N'Z + ?N 2 C - AN c n f 



(23) 



The evolution of the quark mass is determined by eq. (|19|). The solution 
is particularly simple in the MS scheme, where the renormalized quark mass is 
gauge independent. In this case, it can be expressed in the form: 



m MS fa 



.MS 



.MS 



where: 



(24) 



c MS (/*) 



1 / a s fa ' 

2 I 4?r 



(25) 



__\2 _ _ 2 _ _ _ _ _ 

T1-A1T0 + 7 2 + 4i7o - 4i7i - /5 2 7o 



with ft = A/A) and 7, = 7 »/ (2(3 ). 

By using eq. (0), we can convert eq. (|4]) into a relation between the RI quark 
mass at a scale /i and the MS one at a different scale \i: 



c MS (M) Z Rl fa) RI 
c M S fa) ZMS (//0 ) 



(26) 
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This formula can be now easily interpreted by observing that the function: 

„RI 



(A*) 



^ MS (aO c ms 

2 RI (/x) 



(A*) 



(27) 



is the solution of eq. (|19D in the RI scheme. So we can conventionally define a 
renormalization group-invariant mass as: 



ni 



MS 



(A*) 



m 



RI 



(A*o) 



.MS 



(A*) 



,RI 



(A*o) 



(28) 



The mass m is a short-distance quantity, which is both scale and scheme inde- 
pendent. It can be conveniently used in phenomenological applications involving 
quark masses, in alternative to the MS quark mass. 

The relation between rh and the MS quark mass has been computed in ref. |lj]. 
For rif 



3, 4, 5, the result has the form: 



m 



(3) 



m 



(4) 



m 



(5) 



Ot a (ft) ^ 
/ \ -12/25 

a s (fi) 

I \ -12/23 



1 - 



a s (fi) 290 a 2 . /259943 _ 80 

4tt 81 ~ (4tt) 2 V 13122 ~ IP 3 
_ a s (//) 7606 _ of (^) / 446305267 64 
~ 4tt 1875 ~ (4vr) 2 V 21093750 
_ a s (//) 7462 _ (fj) / 344665349 _ 400 
~~ 4vr 1587 ~~ (An) 2 V 15111414 ~ ~23~ 



m MS (//) (29) 







We can use eq. ( |27|) to state our results of eqs. (p~3|)- (P~8|) as a relation between 
the renormalization group-invariant mass m and the mass m . 



m 



(3) 



a s (/i) 



-4/9 



a s (/i) /722 



4tt V 81 
a 2 s (fi) /2521517 536 



+ 2£ 
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(4) 



q a (/i) / 17606 
4tt k 1875 + * 



(5) 



(4tt) 2 V 13122 

/ \ -12/25 

a s (fi) 

a 2 (fi) /3819632767 952 
~ (4tt) 2 V 21093750 15~ 

a s (fi) 



257, 11, 2 

C H C 

81 ^ 6 ^ 



m RI (fj,) 



(30) 



4163 
C3 + 1875 e 
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m RI (/i) 



-12/23 



_ a s (fi) f 15926 
~ 4tt V 1587 
a* (/i) /2559841211 4696 



+ 2^ 
1475 



(4vr) 2 V 15111414 



69 



C3 + 



1587 



m RI (fi) 



RT RT 

We note that in the RI scheme both m and c depend on the gauge, but the 
dependence cancels out in the ratio m. 
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